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^ ! Abstract 
OO , 

■ For a non-Hermitian Hamiltonian H possessing a real spectrum, we introduce a canon- 

■ 

' ical orthonormal basis in which a previously introduced unitary mapping of to a Her- 



mitian Hamiltonian h takes a simple form. We use this basis to construct the observables 
Oa of the quantum mechanics based on H. In particular, we introduce pseudo-Hermitian 
position and momentum operators and a pseudo-Hermitian quantization scheme that re- 
lates the latter to the ordinary classical position and momentum observables. These allow 
us to address the problem of determining the conserved probability density and the un- 
derlying classical system for pseudo-Hermitian and in particular PT-symmetric quantum 
systems. As a concrete example we construct the Hermitian Hamiltonian h, the phys- 
ical observables Oa, the localized states, and the conserved probability density for the 
non-Hermitian PT-symmetric square well. We achieve this by employing an appropriate 
perturbation scheme. For this system, we conduct a comprehensive study of both the 
kinematical and dynamical effects of the non-Hermiticity of the Hamiltonian on various 
physical quantities. In particular, we show that these effects are quantum mechanical in 
nature and diminish in the classical limit. Our results provide an objective assessment 
of the physical aspects of PT-symmetric quantum mechanics and clarify its relationship 
with both the conventional quantum mechanics and the classical mechanics. 
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1 Introduction 

Most of the recent publications on PT-symmetric quantum mechanics focus on the study of 
the spectral properties of yarious (non-Hermitian) PT-symmetric Hamiltonians. The results 
reported in these publications are mainly mathematical. The purpose of the present paper is to 
address some of the most basic problems related to the physical aspects of PT-symmetric and 
more generally pseudo-Hermitian quantum mechanics. In particular, we will offer a complete 
description of the nature and the construction of the physical obseryables and proyide a method 
to compute yarious physical quantities in these theories. We will also elucidate the relationship 
between these theories and the conventional classical and quantum mechanics. 

As our approach is motivated by the mathematical results obtained within the framework 
of the theory of pseudo-Hermitian operators ^1 121 E] ? we begin our discussion by a brief review 
of the relevant developments. 

A central question that arises in connection with the current interest in PT-symmetric 
quantum mechanics [H E| is: "What are the necessary and sufficient conditions for the reality 
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of the spectrum of a linear operator?" Ref. provides the following answer to this question: 
If the operator acts in a Hilbert space Ti and has a complete set of eigenvectors (i.e., it is 
diagonalizable) then its spectrum is real if and only if (one and consequently all of) the following 
equivalent conditions holds. 

(CI) There exists a positive-definite operator^ 77+ : 7i — > H that fulfils 



i.e., H is pseudo-Hermitian JT] and the set^ of all the metric operators r] satisfying H'^ = 
rjHrj'^ includes a positive-definite element. 

(C2) H is Hermitian with respect to some positive-definite inner product (■,■)+ on Ti 
(which is generally different from its defining inner product (-I-).) A specific choice for 
(-,■)+ is {■\r]+-). 

(C3) H may be mapped to a Hermitian Hamiltonian hhy a, similarity transformation, 
i.e., H is quasi-Hermitian j7||H]. 

The framework provided in Refs. ^1121 also explains the connection with PT-symmetry. It turns 
out that, under the same conditions, pseudo-Hermiticity of H is equivalent to the presence of 
an antilinear symmetry, PT-symmetry being the primary example, jSlE]- 

The condition that the Hamiltonian H must have a complete set of eigenvectors may be 
relaxed by extending the analysis of [1113131 to block-diagonalizable linear operators as discussed 
in pm lllj. However, note that physically this condition is intertwined with the requirements 
of the quantum measurement theory. The failure to satisfy it is equivalent to allowing for the 
states that have zero overlap with all the energy eigenstates. As a result, the total probability 
of measuring any energy value for such a state is identically zero, i.e., one can never perform 
an energy measurement on such a state; it must not be possible to prepare it! 

These physical considerations form the basis of a general framework, called pseudo-Hermitian 
quantum mechanics |T2] , that allows for formulating a quantum theory based on an eigenvalue 
problem for a linear operator H acting in a (complex) vector space V . A typical example is 
an eigenvalue (Sturm-Liouville) equation for a differential operator acting in a complex func- 
tion space. Supposing that this eigenvalue problem has a solution, i.e., there are eigenvectors 
ipn G V , one lets Vh be the span of ipn, endows Vh with an arbitrary positive-definite inner 
product, Cauchy completes [13] this inner product space to a Hilbert space 7i, and views H 

"'^An operator is called positive-definite if it is Hermitian and has a strictly positive spectrum. 
^For a discussion of this set, see 
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as a (possibly densely defined) linear operator acting in H.. Then, by construction, if is a 
diagonalizable operator acting in ?i, and the results of [H |21 E] apply. 

As noted in p!l|, the equivalence of the reality of the spectrum of H and the condition (C2) is 
the basic mathematical result underlying the construction of the so-called CPT-inner product 
for PT-symmetric quantum systems Also as shown in JE], one can use the condition (C3) 
to map H to a Hermitian Hamiltonian h acting in Ti. If one identifies the physical Hilbert 
space T^phys of the system with TC endowed with the positive-definite inner product (■,■)+, then 
H and h are unitarily equivalent. 

For models with a finite-dimensional Hilbert space the construction of the Hermitian Hamil- 
tonian h is straightforward. In some cases h has a much simpler form than H, ^H] • The situation 
is quite different for systems with an infinite-dimensional Hilbert space, as almost nothing spe- 
cific is known about the structure of h. It is nevertheless expected to be a generally complicated 
nonlocal (non-differential) operator [T^ . 

The study of systems with an infinite-dimensional Hilbert space is particularly important, 
because it is for such systems that one can seek for an underlying classical system and attempt 
to formulate an associated quantization scheme. Obviously, a proper treatment of these issues 
requires a careful study of the notion of a physical observable in pseudo-Hermitian and, in 
particular, PT-symmetric quantum mechanics. 

It has recently been shown that the formulation of observables in PT-symmetric quan- 
tum mechanics as originally proposed in ^5] and reiterated in JH] is inconsistent with its 
dynamical aspects and that enforcing the rules of the standard measurement theory restricts 
the choice of the observables Oa to linear operators that are Hermitian with respect to the in- 
ner product of the physical Hilbert space Tiphys-'^ Accepting this definition for the observables, 
one can easily show that the unitary mapping that maps H to h also maps the observables 
Oa to the Hermitian operators Oa acting in Ti. This in turn means that a physical system 
described by the Hilbert space T^phys, the PT-symmetric Hamiltonian H, and the observables 
Oa may be equivalently described by the Hilbert space Ti, the Hermitian Hamiltonian h, and 
the observables Oa- 

In this paper we will introduce a canonical basis in which the construction of the Hermitian 

■^To resolve the inconsistency reported in |17| . the authors of |15| have recently revised their definition of 
observables ^]. As noted in [20], it is not known if this corrected definition is generally compatible with the 
requirements of the quantum measurement theory. When the contour defining the boundary conditions of the 
problem is the real line, the definition reduces to ours (and consequently the above compatibility condition 
holds). But even in this case it is a more restrictive definition as it implies that the Hamiltonian must be (not 
only PT-symmetric but also) symmetric, i.e., in x-representation it is a symmetric (infinite) matrix. This leads 
to some undesirable consequences [3D]. 
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Hamiltonian h and the physical observables Oa simphfies considerably. This allows us to de- 
termine the underlying classical system and develop a pseudo-Hermitian quantization scheme. 
We will also introduce and construct the pseudo-Hermitian position operator, the correspond- 
ing position wave functions, and the conserved probability density. As a concrete application 
of our general results we perform a thorough investigation of the PT-symmetric square well 
Hamiltonian, computing the corresponding Hermitian Hamiltonian h, the observables Oa (in 
particular the pseudo-Hermitian position operator), the probability density, the position expec- 
tation values, and the localized states. We will also describe the effects of the non-Hermiticity 
of the Hamiltonian on the latter quantities and discuss the underlying classical Hamiltonian. 

Throughout this paper we will assume that H is diagonalizable and has a nondegenerate, 
real, discrete spectrum. The extension of the results to degenerate spectra may be easily 
achieved following the approach of [H |2l El El- The presence of a continuous part of the 
spectrum does not lead to any serious complications either. For example see j^H 122] • 

2 Canonical Metric Basis 

Let 7i be a Hilbert space and H : Ti ^ Ti he a, diagonalizable linear (Hamiltonian) operator 
having a real, nondegenerate, discrete spectrum. Following 121 El we shall label the 
eigenvalues of H with En and let {iV^n)} denote a basis of Ti consisting of the eigenvectors \ipn) 
ofH, 

H\^^) = E^\^n). (2) 
Then one can construct another basis {|0n)} of H that satisfies jH] 

H^(j)n) = En\4>n), (0n|V^m) = ^mn, ^ |V^n)(0n| = 1- (3) 

n 

In particular, {\ipn), |0n)} form a biorthonormal system and 

H = Y,En\^n){<Pn\, = En\<Pn) {^n\. (4) 

n n 

Here and throughout this paper, for any linear operator acting in ?i, stands for the adjoint 
of A, i.e., the unique linear operator satisfying {■\A-) = {A'' ■ \-). 
A central result of |T] is that the operator 

V+ = 'Y\<Pn){(l)n\ (5) 
n 

satisfies ((H). It is also manifestly positive-definite, because it satisfies 77+ = w^w, where w : = 
Yin \''^){^n\ and {|^)} is an orthonormal basis of Ti, and that it is invertible, with the inverse 
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given by 

n 

We can use ^ to introduce the positive-definite inner product: 

(-,■)+ :=(-|^+-), (7) 

and identify tlie pliysical Hilbert space T^phys witli tlie underlying vector space of Ti endowed 
with this inner product. This means that as complex vector spaces Ti and Tiphys are identical, 
but as Hilbert spaces they are not. 

In view of (PQ), the Hermitian Hamiltonian h of condition (C3) has the form [TB] : 

h = pHp-\ (8) 

where p := is the unique positive(-definite) square root of t]^. The transformation H —>■ h 
corresponds to the linear mapping {tp) —>■ p\ip). It is a simple exercise to check that, for any 
pair l"^), \ip') of state vectors: {ilj,ip')+ = {i^lv+i^') = {p4'\p4'')- Hence as a mapping of T^phys 
onto 7Y, p is a unitary operator."^ 

Now, consider a physical system S that is described by the Hilbert space Tiphys, the Hamil- 
tonian H, and the observables that are Hermitian operators acting in T^phys-^ Because 
p : 7-^phys ^ 7"^ is a unitary transformation, Oq, : 7-^phys '^phys is Hermitian if and only if 
Oq, := pOaP^^ : — > 7^ is Hermitian. This, in particular, means that the observables Oa niay 
be constructed from the Hermitian operators Oq, according to ^Tj 

Oa = p'^OaP. (9) 

Consequently, we can also describe the physical system S using the original Hilbert space 7i, the 
Hermitian Hamiltonian /i, and the observables Oq,. The two descriptions are physically identical 
as there is a one-to-one correspondence between the states and the observables used in these 
descriptions and more importantly the physical quantities such as the transition amplitudes or 
expectation values of the observables do not depend on the choice of the description. 

The main ingredient of the above construction is the operator p = i/r/^. It has three 
important properties: 

PI. As an operator mapping Ttphys to H, it is a unitary operator; 

linear map U : Hi Ti.2 between two inner product (in particular Hilbert) spaces Tii and 7^2 with inner 
products (•, •)! and (•, •)2 is said to be a unitary operator if for all Ci X ^ have (f/(C), U{x))2 ~ {C, x)i7 

. U is unitary if and only if it is invertible (one-to-one and onto) and = [/t. 

^Being a Hermitian operator acting in Tiphys, the Hamiltonian H is also an observable. But as operators 
acting in Ti neither H nor Oa are Hermitian. 
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P2. As an operator mapping 7i to 7i, it is a Hermitian operator; 

P3. As an operator mapping Hphys to Hphys, it is also a Hermitian operator.^ In particular, 
both p and 77+ = are physical observables. 

Property P2 suggests that a natural method for computing the operators h and is to use 
an orthonormal basis {|en)} of 7i that consists of the eigenvectors^ |e„) of rj^. Denoting the 
eigenvalues of r]^ by e„, we have 

T]+\en) = e„|e„), (em|e„) = ^ |en)(e„| = 1. (10) 

These in turn imply 

n n 

In the following we shall refer to {|en)} as a canonical metric basis. 

Let A be a linear operator acting in ?-^, we can uniquely identify A with its matrix repre- 
sentation (Amn) in the basis {|e„)}, where 

Amn •= (Cml^l^n)- (12) 

Because {|e„)} is an orthonormal basis of 7i, the matrix elements of A'^ are given by 

/it — A* (1 Q^ 

mn nm' \ ^/ 

In particular, A is Hermitian with respect to the defining inner product (-I-) oiTi if and only 
if {Amn) is a Hermitian (possibly infinite) matrix, i.e., A*^^ = Anm- 
The following important identities follow from (HJ, ©, © and 

en = Y.\{<Pm\en)?= ^ I (^ni I ^n) , (15) 



m 




hmn \ Hffin • ( ^ ^ ) 



er. 



Furthermore, let o : ^ be a Hermitian operator and O := p op, then 




^This can be easily checked: {■,p-)+ = (•, \r]+p-) = (•|pry+-) = {p ■ \r]+-) = (p-, •)+. 



^Here we suppress the degeneracy labels for the eigenvectors |e„) for simplicity. Note also that in view of 
the non-uniqueness |^ of r]+ one can assume without loss of generality that the eigenvalues of 77+ are 
nondegenerate. 
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Eqs. fjl6p and (fT7j) provide the following expressions for the Hermitian Hamiltonian h -.Ti ^ Ti 
and the observables O : T^phys '^phys- 



h — ^ -i /— -f^mn |em)(en|, (18) 

V ^Tl 

m,n ' 

O = ^ W— |e.^)(e„|. (19) 



m.n 



3 Classical System and Its Pseudo-Hermitian Canonical 
Quantization 

For Ti = L^(R), we can define the rj^-pseudo-Hermitian position (X) and momentum (P) 
operators according to 



-X. . — ^ ^ 4 / Xfnn |Cm)(Cn|; P ■ — ^ ^ ■\ / Pmn |Cm)(Cn|) (^O) 



— X le Ue I P ■=\^ J— 

^"^ m,n V 

where Xmn '■= (em|a;|e„), := {^m\p\^n), and x and p are the usual position and momentum 
operators acting in 7-^ = L^(M). 

Clearly, the T^^-pseudo-Hermitian position and momentum operators satisfy the canonical 
commutation relation 

[X,P]=ihl. (21) 

Indeed, together with the identity operator 1, they provide a unitary irreducible representation 
of the Weyl-Heisenberg algebra which has the physical Hilbert space Tiphys as the representation 
space. The fact that by construction this representation is unitarily equivalent to the standard 
representation of the Weyl-Heisenberg algebra (that has 7i = L^(R) as the representation space) 
is a manifestation of von-Neumann's celebrated uniqueness theorem.^ 

Having introduced the r^^-pseudo-Hermitian position and momentum operators, we can also 
speak of the following rj^-pseudo-Hermitian canonical quantization of classical systems: 



X, p,^P, {.,.}^^-^h-^[.,.], (22) 



where Xc,Pc, and {^-jc stand for classical position, momentum, and Poisson bracket, respec- 
tively. For instance, r^^-pseudo-Hermitian quantization of the classical Hamiltonian for a free 
particle leads to the pseudo-Hermitian quantum Hamiltonian: 

p2 



Hfree = = ^ 

2m ^ 



p2 

2m 



P, (23) 



^This theorem states that up to unitary equivalence the Weyl-Heisenberg algebra has a unique unitary 
irreducible (projective) representation [Oj . 
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which is a generally nonlocal (non-differential) operator. 

Note that in general the Hamiltonian operator H, that is used to construct the metric 
operator and consequently define the above notion of pseudo-Hermitian quantization, does 
not have the standard form P^/(2m) + V{X). For example, a PT-symmetric Hamiltonian of 
the standard form [4^ (with a complex-valued potential f (x)), 

p-\ (24) 

cannot generally be expressed in the form P"^ /2m + V{X) for any real- valued function V. 
Nevertheless, because (in light of property P3) p is also a physical observable, one can express 
p and p~^ and consequently the Hamiltonian ()24|) as certain power series in X and P (modulo 
commutation relations (pT|l .) This in turn implies that the classical Hamiltonian ifc, whose 7]+- 
pseudo-Hermitian quantization yields if, is not generally of the standard (Kinetic-|-Potential) 
type. Rather it is a complicated (non-polynomial) function of Xc and pc- 

The classical Hamiltonian may also be obtained using the Hermitian Hamiltonian h 
which according to (jH)) and (j24j) takes the form 

p'\ (25) 

Again this is a nonlocal operator which can be expressed as a power series in p with x-dependent 
coefficients. This is because (according to property P2) p and p~^ are Hermitian operators 
acting in Ti. The classical Hamiltonian may be obtained by replacing x and p in the expression 
for h by their classical counterparts Xc and pc-, respectively. Clearly the resulting He is identical 
with the one obtained from H. 

Next, we wish to recall a simple procedure for associating a power series in x and p (i.e., a 
pseudo-differential operator) to a nonlocal linear operator K : L^(R) L^(]R). Suppose that 
K may be expressed in terms of its kernel /C : ^ C according to 

{Kij){x) = [ IC{x,x')ij{x')dx'. (26) 

Then for real analytic wave functions, that form a dense subset of L^(R), we can expand tp{x') 
appearing on the right-hand side of in Taylor series about x. Substituting the result in 
ipmi. we find {KiIj){x) = Kijjix) where 

K := 5^(-zn)V(x)^, (27) 

£=0 ^ 
■£ r 

ai{x) = K{x,x'){x' -xYdx'. (28) 



H 



p" 
2m 



+ v{x) = p 



pi 
2m 



+ v{X) 



h = p 



p^ 
2m 



+ v{x) 
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As a result, we have the following (densely defined) identity 

oo 

K = ^a,(x)/. (29) 
If the operator K is Hermitian, we can express (j^^ in a manifestly Hermitian form, namely 

^ oo 

K=-Y,Mx)p' + p'a,{xr]. (30) 

The classical counterpart of this operator is the following real-valued function of the phase 
space (M^). 

oo 

Kc{xc,Pc) ■=^^^[ai{xc)]Pa XcPceM, (31) 

1=0 

where 3? means 'Real part of. 

The results reported in this section clearly generalize to the Hilbert spaces 'HiV) where V 
is or a topologically equivalent subset of M". Together with the results of the preceding 
section, they lead to the following prescription for determining the classical Hamiltonian for a 
pseudo-Hermitian (particularly PT-symmetric) quantum system: 

1. Given the Hamiltonian H, compute a metric operator 77+; 

2. Diagonalize 77+ and construct the corresponding canonical metric basis {|e„)}; 

3. Compute the matrix elements i^mn of H in this basis and use (jl8|) to obtain the Hermitian 
Hamiltonian h\ 

4. Apply the above described method of associating a pseudo-differential operator to the 
operator h, express the latter in a manifestly Hermitian form {h + h^)/2, and take x ^ Xc 
and p Pc in the resulting expression. This yields a classical Hamiltonian He for the 
theory. 

The Hamiltonian He obtained in this way generally involve h. The strictly classical Hamiltonian 
will correspond to evaluating h limit of H^. The latter is an admissible prescription only 
if this limit exists. 

We end this section by making the last step of the above prescription more specific. Using 
ffTHj) . (j^ . and (jnH), identifying the kernel of h with {x\h\x'), and denoting the normalized 
eigenfunctions of 77+ by i.e., 

Snix) = (x|e„), (32) 
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we have 

oo oo 

h = ^a,(x)/, H, = J2^Mxc)]p'c, (33) 

e=o £=0 

aeix) = —y2\—Hmn£rnix) enix'y{x'-xYdx'. (34) 



m,n 



Admittedly, the computation of h as outhned above is too comphcated to be done exactly. In 
Sec. 5, we study its application to a simple PT-symmetric model where a particularly useful 
approximation scheme allows for computing h with any desired accuracy. Finally, we should 
like to add that the above prescription for computing h may also be used to compute the 
pseudo-Hermitian observables such as the position operator X. 

4 Localized States, Position Wave Functions, and the 
Probability Density 

Having introduced the //^-pseudo-Hermitian position operator X we can identify its (generalized 
j25j ) eigenvectors ^^^^ with the localized states of the system. They are defined by 

j^^W^^^W^ VxgR. (35) 
In view of the identity X = p^^xp, we have 

&^ = p-'\x), (36) 
where \x) are the usual position kets satisfying, for all x, x' G M, 

{x\x') = 6{x — x'), / dx\x){x\ = l. (37) 

Jr 

Using these relations and the fact that p : T^phys — ^ is a unitary mapping, we can establish 
the orthonormality and completeness relations for the localized states 

{&\&'^)+ = S{x-x'), [ rfxS(-) = l, (38) 

Jr 

where 'E^^^ denotes the projection operator defined by 

E^-)^ := for all e T^phys- (39) 

Next, consider a particle^ whose state at a fixed time to is described by the state vector 
ip G Tiphys- We can introduce the position wave function: 

^{x):={^^^\i;)+ = {x\p\ij), (40) 



^Here by a particle we mean a quantum system having R as its classical configuration space. 
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and use (|HHj) to expand the state vector if) in the position basis {^*-'^^} according to 

^P= f ^{x)i^''^dx. (41) 

Jr 

As seen from ()40|) . the position wave function is generally different from ip{x). This is a 
direct consequence of the fact that H as an operator acting in 7^ = L^(R) fails to be Hermitian. 
Furthermore, in view of (|HS|) and (jlOl), 

II ^ ||2:= I \m{x)\^dx= I |(^(^),^/>) + |2rfa; = I {iP ^Z^''^) +dx = {^j , i)) + < 



oo. 

Hence as a function mapping M to C, the wave function \1/ belongs to L^(R).^° The converse is 
also true in the sense that every square-integrable function defines a state vector ip G Tiphys- 
Therefore, we may identify L2(R) with the vector space of position wave functions \1/ for the 
system. It is also a straightforward exercise to show that the assignment F of a wave function 
\1/ =: F{ip) to each state vector viewed as a map F : Tiphys -^^(1^), is a unitary operator. 
In order to see this, let V^, G T^phys be arbitrary state vectors and \1/ = F{jp) and $ = F{ip), 
then 

(F(^)|F(0)) = (vI/|$)= / ^ii{xr^x)dx= /(V',e(^))+(e(^),0)+rfx= /(V^,S(^V)+rfa;=(^,0) 



In the following we will assume without loss of generality that ip is normalized with respect 
to the inner product (■, ■) + , i.e., set {ip,4')+ =\\ |p= 1. 

According to the standard quantum measurement theory, the probability of finding the 
particle in a region C R at time to is given by 

nyW:= / \{^'^'\^P)^\'dx. (42) 
Jv 

Hence 

g{x):=\{^^^\ij)^\' = \^{x)\' (43) 

is the probability density of the localization of the particle in space. 

Unlike the naive "probabihty density" g^'^\x) := |-?/'(x)p, g{x) defines a conserved total 
probability. This follows from the fact that H is Hermitian with respect to the inner product 
(■,■)+ of Tiphys- It is instructive to demonstrate the conservation of total probability in the 
position representation. In order to do so, consider the time-evolution of the state vector ijj as 
determined by the Schrodinger equation: 

^hj^m = Hm- (44) 



^"Here we view ^ as an abstract vector belonging to L^(M). The state vector ip also belongs to the Hilbert 
space H which coincides with i^(R). However, these two copies of L^(M) should not be confused. 
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Computing the inner product of both sides of this equation with ^^^^ (using the inner product 
(■,■)+) and employing the completeness relation given in ()38p. we find 

ih^'^{x;t) = h^{x;t), (45) 

where ^(x;t) := (^(^),^(t))+ and h : ^^(R) L'^{R) is defined by 

h<il{x; t) := / /C(x, y)<il{y; t) dy, /C(x, y) := {^^^\ H^^y^)+. (46) 
Jm. 

Because, as an operator acting in Tiphys, H is Hermitian, 

This is sufficient to conclude that /i is a Hermitian operator acting in L^(]R). As a result, 
in the position representation the dynamics is determined by a Hermitian Hamiltonian; the 
time-evolution operator, e~*(*~*o)'^/'^, for the position wave functions is unitary; and the total 
probability 

nR(V'(t)) = / \^ix;t)\^dx 



is conserved. 

The Hamiltonian operator h is directly related to the Hermitian Hamiltonian h. Substituting 
in pUj) and using ((Tj) and (jHl), we have 

}C{x,y) = {x\p~^ri+Hp'^\y) = {x\h\y), 

h^{x;t) = [ {x\h\y)^iy;t)dy= [ {x\h\y){y\pij{t))dy= {x\hpijit)). 



Hence, in light of ()4U|) . the Hamiltonian h is the usual position representation of the Hermitian 
Hamiltonian h, i.e., 

{x\h = h{x\. (47) 

This relationship between the Hamiltonian operators h and h extends to all the physical 
observables. Given an observable O acting in T^phys and the corresponding operator o = pOp~^ 
acting in Ti, we can define an associated Hermitian operator 6 acting L^(M) that realizes the 
action of O on a state vector ip G Tiphys in terms of the corresponding position wave function 
\l/ according to 

0^= [ [6^(a;)]^(^)rfx. (48) 
Jr 

The operator 6 is the position representation of the abstract operator o; 

{x\o = d{x\. (49) 
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In view of (pSj). (pUj) . (|HHj). and (jlOj) . the expectation value of O in a state described by the 
normahzed state vector ip and position wave function \1/ is given by 

/oo 
^(a;)*6*(a;) rfx. (50) 
-oo 

As shown in the preceding paragraphs, one can formulate both the dynamics and the kine- 
matics of the theory using the position wave functions \E'. In this formulation the observables 
and in particular the Hamiltonian are Hermitian operators acting in 7i similarly to the conven- 
tional quantum mechanics. In order to use this formulation, however, one needs a more explicit 
expression for the wave function \E'. We may derive such an expression using the canonical 
metric basis {|e„)}. In view of, dTT)), (El), and we have 

'fix) = y2Uenix), fn:=el/' [ £„(x')*^(a;') rfx'. (51) 

5 Application to the PT- Symmetric Square Well 

The PT-symmetric square well potential, originally introduced by Znojil in [26 , provides a 
simple model with generic properties of general PT-symmetric potentials. Its Hamiltonian is 
given by 



where 



vix] 



= I- + .„(.). (52) 



OO for a; ^ (-f , f ) 

< for xG(-f,0) (53) 



-< for X G (0, 



2 • 

^ 2 ^' 

L G and C ^ Usually one employs units in which L = 2, m = 1/2, and h = 1. This is 
equivalent to using the dimensionless variables 

.-x:=(|)x. p^p:=(|,)p. C-Z:=(|^)C. (54) 
and working with the dimensionless Hamiltonian: 

(55) 



where 







H = 


■ + v(x), 




\ : 


for 


x^(-l,l) 


v(x) = 




for 


XG (-1,0) 




1 -iZ 


for 


XG (0,1). 



(56) 
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In the x-representation, the eigenvalue problem for H takes the form^^ 



+ v(x) - E, 



^n(x) = 0. (57) 



The Hilbert space to which the eigenvectors tpn belong is^^ 

n = L^(M) I ^{x) = for |x| > 1} = G ^^([-1, 1]) | ^^(±1) = O} . (58) 

Clearly, H is not Hermitian with respect to the defining inner product (-I-) of Ti. This is an 
indication that Ti is not the physical Hilbert space T^phys- In order to specify the latter we 
should determine an appropriate metric operator r]^. This in turn requires the solution of the 
eigenvalue equation (j37jl . 

The eigenvalue problem for the PT-symmetric square well admits an essentially explicit 
solution. A detailed discussion is provided in pUllTTj . If Z is below the critical value ^ 4.48 
the Hamiltonian has a real spectrum [^Hl-^^ For these values of the 'non-Hermiticity' parameter 
Z, the Hilbert spaces H and Hphys are identical as complex vector spaces, i.e., they are obtained 
by endowing their common vector space with different inner products. 

In this paper we will only be concerned with the case < Z < Z^. For these values of Z, 
one obtains the following complete set of eigenfunctions of H: 

v.(x) ^ ^'^i-;;' (50) 

[ ?/'„+(x) for X G [0, 1], 

^'^iW := sinh(«:„±) ' ^^0) 

where are arbitrary nonzero real coefficients, 

Kn± = SnT itn, (61) 

:= ^, (62) 
and tn with n G Z"^ are the real solutions of the transcendental equation: 

(Z/tn) sinh {Z/t„) + 2tn sin(2t„) = 0. (63) 



^^The eigenvalues En of the Hamiltonian (|52() are given by 2?i^E„/(toL^). 

^^The Hilbert space associated with the unsealed Hamiltonian H is obtained by changing 1 in (|58|l to L/2. 

^^Z = marks an exceptional point 1^8! where two real eigenvalues cross in such a way that the Hamiltonian 
becomes non-diagonalizable. Once Z exceeds Z-^, H regains its diagonalizability, but a pair of complex-conjugate 
eigenvalues appear in its spectrum [lOL |S] . Increasing the value of Z indefinitely one encounters an infinite 
number of exceptional points passing each of which produces a complex-conjugate pair of eigenvalues. 
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The eigenvalues E„ are given by 

E^ = -{Kl^ + iZ)=tl-sl. (64) 

Usually the coefficients a„ are fixed arbitrarily [201 or kept as unimportant free coefficients 
|27j . We will fix them in such a way that in the limit Z — ^ 0, the eigenfunctions ipn of ()59p 
tend to the well-known normalized eigenfunctions ijjn^ of the conventional (Hermitian) square 
well Hamiltonian (the case Z = 0): 

hmV'„ = V^W:= V'nl,^^. (65) 

Because, by construction, ipn are also the eigenfunctions of the PT operator, the continuity 
requirement (j65j) constrains ipn^ to be PT- invariant. The normalized and PT-invariant eigen- 
functions of the Hermitian infinite square well potential {Z = 0) are, up to a sign, given by 

^(0) (x) = z^" sin (x + 1)] , (66) 

where 

1 + (-1)" 

f^n := ^\ ^ . (67) 

The eigenfunctions ()66p form an orthonormal basis of the Hilbert space ()58p. We will denote 
the corresponding abstract basis vectors by \n), i.e., 

(x|n) := ^^°)(x), \/neZ+. (68) 

The continuity requirement (j65j) together with Eq. (j66|) restrict the coefficients a„ of the 
eigenfunctions (jSHI)- Specifically, if we only keep the leading order term in powers of Z and 
neglect the higher order terms, we find 

(-l)LtJ(A) , (69) 



where L|J stands for the integer part of |. Both the eigenvalues E„ and the eigenfunctions ipn 
are therefore determined once we obtain the solutions of (jUHj) . As shown in [21], this equation 
may be easily solved numerically for various values of Z. In this paper, we will solve this 
equation perturbatively by expanding the relevant quantities in powers of Z. 

5.1 Perturbative Calculation of i/jn and En 

Suppose that t„ admits a power series expansion about Z = 0: 

oo 

^n = $^e^^^ (70) 

k=0 
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where tn'' G M are to be determined. Substituting (fTOj) in (jU^ . expanding both sides of the 
resulting equation in powers of Z, solving it term by term for tn \ and using ((70)), (jH^ and 
(jMl) . we find 



En 



/Txn\ 



(-I)"7r2n2 



3 + 



6 



1 + 



-1) + 



4 + 
1)"] . 



(-I)"7r2n 
6 



:-l)'^7rV 



where 



5 + 2(-l)^ 



2Z 



V :- 



{liny 



(71) 
(72) 
(73) 

(74) 



Tcn 



,2 14 



and (9(i/'^) stands for terms of order and higher. 

Eqs. ()71|) - (fTSj) reveal the curious fact that the effective perturbation parameter is 2Z/ ( 
This is a clear indication that the non-Hermiticity of the Hamiltonian H only affects the low 
lying energy levels. This property of the PT-symmetric square well Hamiltonian — which has 
been previously known [22] — is particularly significant, for as we explain below it implies that 
within the confines of the perturbation theory all the infinite sums appearing in the expressions 
(0), (fTH|) . and ()5ip for the metric operator rj^, the Hermitian Hamiltonian h, and the position 
wave functions may be safely truncated. For example, for Z = 1, u ^ 0.2/n^. Therefore, if we 



set Ei°^ 



E, 



n 1 1^=0 



(7m/2)^ and use q„ 



iFjnfEn^ — 1| as a measure of the contribution of the 
non-Hermiticity of the Hamiltonian to the energy eigenvalues, we find for n > 10: u < 2 x 10^^ 
and g„ < 1.3 X 10"^.^^ More generally, we can ignore the effects of the non-Hermiticity param- 
eter Z for all the computations involving the levels with n > 10 and still obtain results that 
are accurate at least up to three decimal places. 

In the following, we will employ an approximation scheme that neglects the effects of the 
non-Hermiticity parameter Z for all levels with n greater than a given number A^. In view of 
the above discussion, the results obtained using this approximation will have an accuracy of 
the order of 

We will respectively refer to and as the 'order' and the 'accuracy index' of our approxi- 
mation scheme. 



1, 



^''The condition that the above perturbative calculations would be unreliable for the ground state, i.e., v •■ 

corresponds to Z w 4.92 which is slightly above the critical value Zi, = 4.48. 

-"^^The value Z = 1, say for an electron (m « lO"'^" Kg) confined in a nanometer size well [L w lO^^m), 

corresponds to an energy scale C ~ 0.1 ev for the potential H53() . This is comparable with the ground state 

energy [Ei « 0.5 ev) of the corresponding Hermitian infinite square well potential. 
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5.2 Construction of a Canonical Metric Basis 

Having obtained the expression for (f7T|) and (f?^ for t„ and s„, we can compute and use 
Eqs. (j^IJj) . (jUUj) . and (I^Uj) to determine the eigenfunctions ^/'n of the Hamiltonian H as a power 
series in v. 

The computation of a metric operator r^^, however, involves the eigenfunctions (pn of the 
adjoint of H. It is easy to see that H^^ = H|^^_^. This suggests that 

Xn ■■= i^n\z^-z (76) 

are eigenfunctions of H"''. The eigenfunctions that together with %pn form a biorthonormal 
system for the Hilbert space 7i, are obtained by properly normalizing Xn- They are given by 

0„(x) = iV-iXn(x), (77) 

where 

Nn ■■= i'ipnlXn) = / ^„(x)*X„(x) 



-1 

tn sin(2t,i ) [cos (2tn ) — cosh(2an )] 



1 — cos(2t„) cosh(2s„) + 



[cos(2t„) — cosh(2sj 



(7J 



and 



Next, we construct the metric operator (0) using the approximation scheme described in 
the preceding section and the orthonormal basis {\n)} consisting of the eigenvectors (j68p of the 
ordinary Hermitian infinite square well. 

In the iV-th order approximation, we have 

N N 

liP^) ^ ^ |n), for all n > iV, (80) 

N 

where ' ~ ' stands for an equality that is valid up to terms of order z/^v- Combining (j8(J|) with 
(jS)) we have 

N oo 

E \n){n\ = l + 6vi''\ (81) 



N 



n=l n=N+l 
N 



5r/f^ := ^(|0.)(0.|-|n)(n|). (82) 



n=l 



This relation shows that the metric operator 77+ is essentially determined by its projection onto 
the span Ti^v of \n) with n < N. Indeed, at this order of approximation, Ti^ may also be 
identified with the span of ipn with n < N, 01 the span of 0„ with n < N. 
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Clearly, the Hilbert space Ti is the direct sum ofHN and its orthogonal complement Hj^ 
{( G 7^1 (CIV') = OjWip G Hn}- As shown by (|HT|) both of these are invariant subspaces 
of ?7+. Therefore, one can solve the eigenvalue problem for 77+ by restricting it onto and 
Hj^ and diagonalize the resulting operators separately. The restriction of 77+ onto Hj^ coincides 
with that of the identity operator. In particular, it is diagonalized in the basis {\n) : n > N}, 
and we have 

N N 

|e„) ^ \n) and e„ «i 1, for all n > N. (83) 

The restriction of rj^ onto TYat yields a Hermitian operator having a Hermitian matrix 
representation S in the basis {\n) : n < N}. According to (|HT|). the matrix elements of S are 
given by 

N 

£mn = {m\r]+\n) = ^(m|0fc) (0fc|ra), (84) 

k=l 

which in view of ^ - ^ - jlll), and {721) can be computed explicitly. The 

computation of the metric basis vectors |e„) with n < is equivalent to the diagonalization of 
the Hermitian matrix S. The latter can be done both numerically and perturbatively. 

Let {e„} be a set of orthonormal eigenvectors of £ so that Sen = ^nCn- Then clearly, up to 
permutations of the labels, e„ coincide with e„ for n < N. The canonical metric basis vectors 
|e„), with n < N, are also related to the eigenvectors e„. To make this relation explicit, we 
introduce the unitary N x N matrix U whose columns coincide with the vectors e„. Then, in 
view of (fTTj) . it is a straightforward exercise to show that 



N 

|e„) ~ Umn\^) and e„ ^ e„, for all n < N. (85) 

m=l 

Clearly, the above approximation scheme would be consistent only if in the above calculation 
of Smn, and e„ one takes into account the contribution of the terms of order for which 
T-'^ > ^N- We can use Eqs. (|74j) and (j75|) to make this condition more explicit. A simple 
calculation shows that the negligible terms are those of order with £ > £„, where 

^ ^ - - In ^= ^ 0.789 - — -. (86) 



''' 1 1 " 1 / — — I 

mra + r \y2Z J 

Clearly, for smaller values of n one should include higher order corrections, the highest order 
term being of the order v^^ . For example, for Z <1 and < 25, ii < 5.03 and one can safely 
ignore (9(z/®). For the values Z = 1 and = 25, the results will have a minimum accuracy of 
the order of z/jv ~ 3.2 x 10~^. Similarly for Z < 0.5 and < 100, £1 < 5.02 and one can again 
ignore (9(z/®). The minimum accuracy corresponding to the values Z = 0.5 and A^ = 100 is of 
the order of un ~ 1.0 x 10~^. Note that for all the above values of Z and A^ the terms given 
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explicitly in Eqs. ()71|) - ()73p are sufficient to perform a consistent perturbative calculation. A 
direct check of the validity of this statement is to compare the exact value E^^'^'^* of E„ obtained 
by an accurate numerical solution of (jUHj) and the perturbative value E^^''*' of E„ calculated 
using ()73p and ignoring 0{i>^). Clearly, the largest difference is for n = 1. If we express 

j^exact = ^ (1 + e^^*), Ef *■ = Y (1 + er*0, 

we find, for Z = 1, ef^^^ = 0.0415652, e^^'^^' = 0.0415527. This is in complete agreement with 
our expectations, because the difference, e™'^* — e^^^' = 1.24465 x 10~^, is much smaller that 
the accuracy index i'n=25 ~ 3.2 x 10""^. 

Next, observe that our approximation, in particular (jSlj) . corresponds to ~ r^^'' where 

\ ^ if 7/^ e nj,. 

It is a reliable approximation, only if the distance between t]^'' and the identity operator 1, as 
defined by 



aN := Y trace [(r^Y^ - 1)2] = trace [(^r^y'O'], (87) 

has a finite large A^-limit. This makes a useful measure of the validity of the above ap- 
proximate calculation of the canonical metric basis {|en)}- Clearly, is just the Frobenius or 
Euclidean distance jSO] between S and the N x N identity matrix J: 



(Tat 



£-1 \\2= ^/tiace[{£ - If 



■ N 

n=l 



n 1/2 



(88) 



where || ■ ||2 stands for the Frobenius (Euclidean) norm jSU], defined for every square matrix 
M by 

II M ||2:= Vtrace(MtM). (89) 

Fig. ^ shows a plot of for Z = 1 and < 25 for which we can neglect 0{v^) in our 
calculations while retaining the consistency of our approximation scheme. The graph of (Tn 
clearly shows the desired behavior even for 20. The effective slope: ctn — cat-i of the 

graph has the values 4.4 x 10~^ and 2.5 x 10~^ for N=20 and N=25, respectively. 

Next, we describe another way of checking the reliability of our approximation scheme. In 
view of (jHUj) and the matrix elements 

:= (m|H|n) (90) 



may be approximated as 



HW^^ -'-^^ (91) 

J TT n Omn lOr 771,71 > JM, 

20 
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where 



Figure 1: Graph of (Tn as a function of N , for Z = 1. 



N 



■.= Y.^k{m\i„){ct>k\n) 



(92) 



fc=i 



are computed by substituting (fTTj) and (f7^ in (jHH) and using (fTSj). and (f77j) . 

Consider the N x N matrix X^'^^ with entries (jU^ and let J^*-^-* be the N x N matrix with 
entries 3^,^^ := Hmn for all m,n < N. Then 3^*^^) may be computed exactly using (j^ . 



and (|9(jp. whereas the calculation of X^^^ uses our approximation scheme. In order to 
compare A'(^) and 3^(^\ we will first introduce the normalized matrices: 



y 



and use the Euclidean distance between the matrices X^^"^ and y^'^\ namely 



(93) 



as a measure of the accuracy of our approximation. A reliable approximation corresponds to a 



vanishing large- limit of Sjy-^^ Fig. [21 shows the plot of S^r for < 25 and Z = 1. Tabled 
lists the values of S^r for various values of A^ and the corresponding values for the accuracy 
index u^. The results indicate that even for A^ = 10 we have a highly reliable approximation. 



"'^^Our use of the normalized matrices ^^"^^^ and CP^^-* stem from the fact that H is not a bounded operator. 
It allows for the interpretation of the term 'a vanishing large- limit' as 'S^v «C 1 for sufficiently large N\ In 
practice this means Sat < vn- 
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Figure 2: Graph of Sat of Eq. as a function of A^, for Z 
Note that S25 ^ 1.1 x 10^^. 



1. 



N 


Sat 


Sat 




10 


8.2 X 10-5 


1.7 X 10-5 


2.0 X 10-3 


15 


3.9 X 10-5 


9.0 X 10-6 


9.0 X 10-^ 


20 


1.7 X 10-5 


3.4 X 10-6 


5.1 X 10-^ 


25 


1.1 X 10-5 


2.4 X 10-6 


3.2 X 10-^ 



Table 1: Values of Sat of Eq. (jHSl), Sn of Eq. (jHHj), and the accuracy 
index Vf^ = 2Z/ (tcN)"^ for Z = 1 and various relevant values of N. 

5.3 Construction of the Hermitian Hamiltonian h 

We can use the above approximation scheme to compute the Hermitian Hamiltonian h := pHp-^ 
(respectively h = pHp~^ = 2^^h/(mL^) ) that is associated with the PT-symmetric square well 
Hamiltonian H (respectively H). In order to do this we first use (fTHj) . (jHUj) . (jH^j) . (jH3j) . to express 
h in the form 

N I 00 

h~ ^ /^H„„|e^)(e„|+ Yl hW|^)(^I + 

m,n=l ' " m,n=N+l 

= H + (5H, (94) 



™ — 1 — AT I 1 
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where Hmn is given by 

[ (e^|H|n) ^ (e™|H|V^„) = E„(e^|n) ^0 Vm < A^, n > A^, 

^ / / — \ 
'■= E J— amn|e^^)(e„|-H(^),|m)(n| . (96) 

m,n=l ^ ' ' 

Substituting (jHSI) in (jHEI) and using ()85|1 and 



we find 



m=l 



AT 



S}i= J2 SR^lHinl (98) 



m,ri=l 

where, for all m, n G {1, 2, ■ ■ ■ , A^}, 



AT 



^^mn ■~ ^mj' ^jj^,^ ^ H^^. (99) 

i,fc=l 

To confirm the consistency of our approximate calculation of h, we check its Hermiticity. 
To do this we compare the N x N matrices Q^^^ and Q^^^"!" defined in terms of their entries 
according to Qmn '■= {m\h\n) and Qmn'^ '■= Qnm = (n|h|m)*. Clearly the condition h"!" = h is 
equivalent to 

QiN)i_^QiN) N^OO. (100) 

Noting that H is not a bounded operator, we follow the method of the preceding section 
and define the normalized matrix Q^^^ := Q^^^ / \\ Q^^^ ||2. This allows us to identify the 
Hermiticity condition (|100|) with 

Sn :=|| - Q^^)^ II2 < i^N. (101) 

Fig. El shows the plot of Sn for < 25 and Z = 1. Table [T] shows some typical values of Sn- 
The results depicted in Fig. Eland Table [Hare in complete agreement with (jl01|) . 

Having obtained the matrix elements of SB. in the basis {|^)}, we can compute its integral 
kernel, 

/C,Jx,x'):=(x|5H|x'), (102) 

and the corresponding pseudo-differential operator in the x-representation. In view of 
(jUH|) and (pg|). 

N N 



}CsnM= E '^H(^^V^^°)(x)V^f (x')*= E ^-r^sm ^(x + 1) Sin ^(x' + l) , (103) 



m,n=l m,n=l 



2 



nn 



2 
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Figure 3: Graph of S^- of Eq. (jlUip as a function of A^, for Z = 1. 
Note that S25 = 2.4 x 10"^ 



where 



(104) 



Next, we follow the derivation of Eqs. ()27|1 - (j29jl to express SH as a series in powers of the 
momentum operator p. This yields 



6ii = Y,H^)p' 



(105) 



i=0 



where 



1 



(106) 



N 



- I /C,jj(x,x')(x' -x) 

1 



■dx' = 5^ A, 



sm 



m,n=l 



vrm 



— (x+1) P„,(x), (107) 



sm 



-KTl 

lT 



(x' + l) (x'-x)W. 



Using ^ - dSni) and ^ we then obtain 

00 00 
h ~ H + 5^ 5,(x) p^ = p' + v(x) + W P' 



(109) 



The analogous expression for the Hermitian Hamiltonian h associated with the unsealed Hamil- 
tonian H is 



^ + v{x) + 7^(x) / 
£=0 



where 



7^(x) := 



L'-%{2x/L) 



(110) 



(111) 
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The integral in p08|) may be evaluated analytically. A simple change of variable reduces 
it to an integral of the form y"^ sm{ny)dy that may be looked up in |31j. Substituting the 
value of this integral in (jlOSj) and doing the necessary algebra, we find 

e 

P„,(x) = 5^an,Kx+l)'~^ (112) 

j=0 



where 



_lY2j+l \^ (^_iy+n+k+l (_l)L§J[l + (_l)i] 



(113) 



and [|J denotes the integer part of |. As seen from fjll2p . Vni is a polynomial of degree i. 

Using (|112j) we can obtain a more explicit expression for the coefficient functions 6i appearing 
in (ITnn|l and (ITTT|l . Substituting (|TT^ in (imTj) and introducing 



TV 



n=l 

we have 

e N 



Seix) = ^ ^ brake sin 



— (x + l)J(x + l)^-^ (115) 

k=0 m=l 

This relation together with (jlOfij) yield the desired expression for ^^(x). 

Next, we recall that the standard calculation ()50|) of the energy expectation value for a state 
vector ip G T^phys, which uses the position wave function \l/ introduced in Sec. 4, involves the 
representation h of the Hermitian Hamiltonian h: 

H^)+ = (ij, H^)+ = ^ J ^ ^ hvl/(x) dx. 
Using (x|h = h(x|, the identity 

/ifln = — i 

dx 



(x|p = -z-^(x|, (116) 



and Eqs. ^M^ . (^11^ . (fTT^ and (fTTIjl . we have 



6S-5j, + vW + 5:mx)5j,. (117) 

£=0 

Figs. Eland El show the plots of the real and imaginary parts of (5^(x) for £ = 0, 1, 2, 3, Z = 1 and 
N = 20. As seen from Fig. the graph of '^{So) is reminiscent of the approximation of a step 
function, namely iv, with the first few terms in its Fourier series expansion. In the appendix 
we offer an explanation for this observation. Furthermore, these figures suggest that Se{x) and 
consequently ^^(x) have a vanishing large-£ limit. This can be established analytically. Using 
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ffTUHl) - (iniHl) and the fact that for all x,x' e [-1, 1] both (x - x')/2 and sin[7rm(x' + l)/2] are 
bounded by 1, we can easily conclude that 

M2^ 

I^.WI < (118) 
where M is a positive number depending on A^n- This shows that, for all x G [—1,1], 



5.4 The Classical Hamiltonian 

Having obtained the Hermitian Hamiltonian h for the PT-symmetric square well, we can use 
the prescription described in Sec. 3 to obtain the following expression for an underlying classical 
Hamiltonian. 

H,{xc,Pc) ^ ^ + J2^be{xc)]pi (119) 

e=o 

where Xc € [— L/2, L/2] and pc € M. ()119p is a meaningful relation, only if the series on its 
right-hand side converges. In view of ()106|) and (jlllll . the latter is equivalent to the convergence 
of 

oo 

J2Wi{^cM, (120) 

e-0 

where x^ G [-1, 1] and Pe G R. According to (ITTH|1 . for all x^ G [-1, 1], |SR[(5£(xJ]| < M2^/i\. 
Using this relation (and performing the comparison and ratio tests [32j) we can easily show 
that the series p20|) converges (absolutely) for all values of Pc G M and x^ G [—1,1]. Hence 
(|119|) is a meaningful expression yielding a well-defined classical Hamiltonian. 

Notice that He as given by (jll9|) is not the classical Hamiltonian in the strict sense that it 
would not involve h. This is simply because we have not evaluated the h ^ limit. Indeed, 
for the PT-symmetric square well Hamiltonian the assumption that this limit exists has 
drastic implications. This is simply because, according to (j33)), if we assume that the coupling 
constant ( appearing in does not depend on h, then taking the limit h ^ corresponds 
to Z ^ oo. This implies the occurrence of an infinite number of complex eigenvalues which in 
turn indicates that the system does not admit a unitary quantum mechanical description 6 . 
The only way in which one can retain such a description and at the same time be allowed to 
take the limit — > is to assume that ( depends on h and is at least of order h"^. 

We can reach the same conclusion by noting that the condition Z < Zj, ^ 4.48, for the 
possibility of formulating a unitary quantum theory for PT-symmetric square well, is equivalent 
to C < 8.96/i^/(mL^) < 2Ei < E2, where Ei and E2 are respectively the ground state and first 
excited states of the system. For a light molecule, say O2, with mass m ^ 30 GeV, confined in 
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Figure 5: Plot of Q[Se\ for i = 0,1,2,3, Z = 1 and N = 20. The 
dashed curve is the graph of iv. The resemblance of the graphs of 
'^[Sq] and iv is described in the Appendix. The difference between 
Q[S£] with £ > 1 are too small to be distinguished in the energy 
scale determined by the potential which is unity in the units used. 
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a micron size well - which should allow for a classical description - we find ( < 1.3 x 10~^^ eV. 
This corresponds to the classical molecule moving with a speed v < 4.7 x 10^^*^ m/s and a 
temperature of T < 1.6 x lO"^'* K. These numbers provide a conclusive evidence that non- 
Hermiticity effects quantified with the coupling constant C are quantum mechanical in nature 
and have no classical counterpart. 

Furthermore, recall that the classical limit /i — is meaningful if it is accompanied with 
taking n — oo in such a way that hn stays constant.^'' But as we explained in subsection 5.1, 
in the limit n — > oo the effects of the non-Hermit icity of the Hamiltonian disappear. Therefore, 
the classical limit of all the theories with different allowed values of Z < Z^, coincides with that 
of the Hermitian infinite square well {Z = 0). 

The above discussion of the classical limit of the PT-symmetric square well Hamiltonian is 
based on the requirement that the corresponding quantum theory has a h —>■ limit. This is the 
conventional way of defining the classical limit of a quantum system. Yet we can consider the 
/i-dependent classical observable He and view it as a classical Hamiltonian with the property 
that its pseudo-Hermitian quantization with appropriate (symmetric) factor ordering yields the 
Hermitian Hamiltonian h. 

There is also another approach for determining a classical Hamiltonian for PT-symmetric 
quantum systems EH] • It involves a direct replacement of the operators x and p, that appear 
in the expression for the quantum Hamiltonian operator H, by the classical position Xc and 
momentum pc and letting the latter take complex values. If one applies this prescription to the 
PT-symmetric square well and enforces the condition of the existence of a proper ^ — limit, 
then again the condition Z < Z^, implies C — ^ 0, and one recovers the classical Hamiltonian for 
a free particle confined in an infinite (real) square well. However, if one does not identity the 
substitution p ^ Pc and x ^ Xc with taking h ^ in (|^. then one obtains a complex- valued 
'classical Hamiltonian', namely 



It is the classical Hamiltonian dynamical systems defined by such complex 'classical Hamil- 

tonians' that are studied in O EHl- Although we acknowledge the interesting mathematical 

consequences of this study and its relevance to the use of complex WKB approximation in 

calculating the energy levels of various PT-symmetric models |34j, we are inclined to adopt 

the standard definition of a classical observable which requires the latter to be real-valued 

|36j.^^ The PT-symmetric quantum mechanics also makes an implicit use of this definition in 

-'^''This follows from H73(l and the requirement that in the classical limit not all the energy levels collapse to 
zero. 

^^This is because we are not aware of any other precise definition of a classical observable. 
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insisting that the eigenvalues of the observables, in particular the Hamiltonian, be real . 
According to this definition, the observables Xc and Pc assume real values, and H'^, which is 
a complex-valued function of Xc and pc, is not a physical observable. In particular, it cannot 
serve as a physical classical Hamiltonian (for a system with a one-dimensional configuration 
space) . 



5.5 Construction of the Observables 

The construction of the observables O : Hphys — ^ '^^phys for the PT-symmetric square well 
mimics that of the Hermitian Hamiltonian h. We begin our calculation of O by employing our 
approximation scheme to express ()19|) in the form 

N ,— 



m,n=l * m,n=N+l 



N oo 



){n\ + ye^o„m|n)(e^ 



m=l n=N+l 

where o : 7i ^ 7i is a Hermitian operator. 



N 



+ So, (122) 



^f,k=i^Ljofk^kn for m,n<A^, 

^ Ef=iWl,oS for m < AT, n > AT, 



{en.\oK)^i , ^^■nj:V\,rfn. . (123) 
o^^l := (m|o|n), (124) 

Af AT oo 

6o := ^Am„|m)(n| + ^ ^ {Bmn \m) {n\ + Cnm \n) {m\) , (125) 

m,n=l m=l n=A''+l 

TV 

^mn ■— / ^ <^m] ^jk '^kn '-'mn^ K^"^^) 
j,k=l 
N 

Bmn := (127) 

k=l 
N 

Cnm ■= y^ojfc i^km -^km), (128) 
k=l 

Smk stands for the Kronecker delta function, and we have used (jHHj) . (jH3j) . and (j^ . 

We can express 60 and consequently O as power series in p with x-dependent coefficients 
similarly to our derivation of piOj) . Note however that in this case we have to deal with the 
infinite sum appearing in ()125|) . The presence of this sum is a manifestation of the fact that 
(unlike the Hamiltonian) a general observable will mix the state vectors \n) with n < N with 
those with n > N. 
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The power series expansion of O in powers of p has the form: 

oo 

O ^ o + ^cu,(x)p^ uj,{x) =t'Coeix), (129) 

I oo 



A:=0 m=l 



2]n=l^'''" '""C'mnOnfc^ for m > N, 

(132) 

where a„fc£ are the coefficients given in (|113p . Again in calculating the expectation value of O 
for a normalized state vector G Tiphys with position wave function \E', we use the position 
representation 6 of o (defined by (x|o = 6(x|): 

{ip,Oilj)+ = j ^(x)* 6^(x) c/x. (133) 

Note also that by replacing (o, x, p) in (jl29|) with the corresponding classical quantities (oc, Xc, Pc), 
we obtain a generally complex- valued function fic(xc,Pc) (provided that the corresponding in- 
finite series appearing in ()130|) and ()131|) converge.) Clearly, taking o = h, we have flc = H'^, 
where H'^ is the complex Hamiltonian ()121|) . 

In order to compare the operators o and O we represent them in the ordinary position 
representation, i.e., compare 6 with O. The latter is defined by (x|0 = 0(x|. Using ()130p . 
(HnH), and (ITTH|) . we have 

oo ,1 

d = d + J2^d-)^r (134) 

£=0 

A concrete example is the dimensionless position operator X := 2X/L: Letting o = x, we 
find 

oo 

X = x + J]cuf^(x)p^ (135) 

e=o 

where uj^^\x) denote the value of uji{x) obtained by setting o = x in (jl24|) . The x-representation 
of X has the form 

X = x + ^i;f)(x)— , (136) 
where ujf^\x) is the value of ue(x) for o = x. Notice that the infinite series in (11311) that defines 



a)^(x) converges quite rapidly. This allows us to obtain an approximate value for this series (for 
any value of i) by summing just the first few terms. We can include enough terms in this series 
so that the approximation error becomes smaller than our accuracy index z/^r. Figs. IHl and [7| 
illustrate the plots of the real and imaginary parts of uJi^\ for i = 0,1,2,3, Z = 1, and = 20, 
that we have obtained in this way. 
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Figure 6: Graph of ^Cof^] for £ = 0, 1, 2, 3, Z = 1, and = 20. 
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Figure 7: Graph of '^\^f\ for £ = 0, 1, 2, 3, Z = 1, and iV = 20. 

5.6 Probability Density, Position Measurements, and Localized States 

According to PHj) . the probabihty density for the locahzation in space is given by the modulus 
square of the position wave function We can employ our approximation scheme to reduce 
the expansion (jSlj) of into the finite sum: 



N 



N \ — V 

\E'(x) ^ ipi^) + / , ctn sin 



n=l 



TTfl , 



(137) 



where 



N 



i^^ f^l'^ /"(o) - /■(' 

'' I / J ^nm J m J n 



(0) 



\m=l 



1 
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vrm 



L 2 



(138) 
(139) 



and we have made use of ^ - (jHEl), (|HS1), (|HS1), (EH), and (j^Tjl . 

Having obtained the general expression for the position wave function, we can compute 
the probabihty density g{x) := |\l/(x)p. Note, however, that the latter expression is valid for 
the normalized wave functions. Fig. |H1 shows the plots of the difference Ag of the probability 
density g for ^(x) = A/; sin[n7r(x + l)/2], with n = 1,2, for = 20 and Z = 0.3,0.7,1 
with that (qq) for Z = 0. Fig. El gives the plots of the probability density difference Ag for 
V^(x) = Mn sin[n7r(x + 1) /2] with n = 3, 4, • • ■ , 8, = 20, and Z = 1.^^ 

We can also use the expression p37|) for the position wave function to compute the position 
expectation value: 

{ip,XiP)+ = j x|^(x)|2dx. (140) 

and the uncertainty in position 

Ax = y/(^,X2^)+-(7A,X7A)^ (141) 

TableHgives Ax for V^(x) = sin[n7r(x+ l)/2] with n = 1, 2, ■ ■ ■ , 7, Z = 0, 0.5, 1 and = 20. 
It turns out the calculation of the same quantities using = 10 yields results that differ from 
those listed in Table El by numbers that are smaller than 10~^. This is another confirmation of 
the consistency of our approximation scheme. Furthermore, note that as we expect the effect 
of the non-Hermiticity of the initial Hamiltonian (j52j) diminishes as n increases. Already for 
n = 6, its contribution to position uncertainty is smaller than the accuracy index i/^o = 0.0005. 

As seen from fjl35p . in the x- representation the position operator X is a pseudo-differential 
operator. This in particular means that the expectation value of X in a state described by 
the state vector ip depends on all the derivatives oi ip. In this sense unlike the usual position 
operator, X is a nonlocal operator. Note however that this nonlocal character of X manifests 
itself only if one insists on using the usual position representation il){x) of the state vectors ip. 
This is not a reasonable choice, because being a non-Hermitian operator acting in Tiphys the 
usual position operator x is not a physical observable. 

Probably the best demonstration of the nonlocal nature of X is provided by the shape of 
the position state vector that is localized at ?/ G {—L, L). As a function belonging to it 
has the form 

oo 

i^y\x) = {x\i^y^) = {x\p-'\y) = Y,^n'^'e4x)e4yr, 

n=l 

where we have employed (jHUj) . (fTTj) . and (jSH). Using the same method as the one leading to 
^^Here A/"„ are appropriate normalization constants. 
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Figure 8: Graph of — Q — Qo for '0(x) = jV„ sin[n7r(x -|- l)/2] with n — 1,2, 
N — 20 and Z — 0.3, 0.7, 1, where g is the probabihty density and is its value for 
Z — and are normahzation constants. 
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Axl 7—1 — Axl 7— n 
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n nni 
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u. oouo 






o 




5571 


555Q 


—0 001 6 


4 


0.5663 


0.5660 


0.5652 


-0.0011 


5 


0.5704 


0.5702 


0.5697 


-0.0006 


6 


0.5725 


0.5724 


0.5721 


-0.0004 


7 


0.5738 


0.5737 


0.5735 


-0.0003 



Table 2: The position uncertainty Ax for ip{x) = NnSm[mT{x + l)/2] with 
n = 1, 2, ■ ■ ■ , 7, Z = 0, 0.5, 1, and = 20. Note that for n = 6 and 7 the 
difference between values of Ax for Z = 1 and Z = is smaller than the 
accuracy index z/20 = 0.0005. 

(jl37j) we can express this relation as 
where 



N 



-1/2 

mn 



5n 



nm , L. 
[x + - 



sm 



Tin L, 



(142) 



m,n=l 

and we have used (jHZI)- 

Fig. ^1 shows the real and imaginary parts of for L = 2, AT = 20, Z = 1, and 

y = —1/2,0,1/2. Fig. ITD shows the real and imaginary parts of ^^^^^\x) for L = 2, N = 20, 
and various values of Z. As expected the spreading of the localized state is an increasing 
function of the non-Hermiticity parameter Z. 



5.7 Dynamical Consequences of Non-Hermiticity 

In the preceding subsection we discussed the computation of the observables and the associated 
physical quantities. These provide information on the kinematical content of the PT-symmetric 
square well. In this subsection, we investigate its dynamical content. 
The time evolution of an initial state vector ipito) G Hphys is given by 

^(t) = e-'^'-''^"'^ij{to), yt e R. (143) 

Alternatively, in terms of the dimensionless time parameter: 




(144) 
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Figure 10: Graph of ^^^y^] and ^[C^^)] for L = 2, = 20, Z = 1, and y = 
— 1/2, 0, 1/2. Note that 3ft [^*^^^] has a 5-function singularity at y and that except for 
this singularity the scale of variations of is much greater than that of 3ft[C^^'']- 
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we have 

V^(r) =e-*(^-^°)»^(ro), (145) 

where tq = 2hto/{mL'^). 

Expanding iP{t) in the basis {ipn}, we can express fll45|) in the form 

oo 

ij{T) = 5^c„e-^(^-^«)'^"^/.„, c„ := {(Pnmro)). (146) 

n=l 

Next, we employ our A^-th order approximation scheme. Using ()8U|) . we have 

N oo 
n=l n=Af+l 

where 

:= (^i°)|^(ro)), (148) 

and we have used E„ ^ E^'''' := n'^n'^/A for n > N. 

To explore the dynamical effects of the non-Hermiticity of the Hamiltonian (for Z 7^ 0) we 
compute the position expectation value for the evolving state vectors iP{t) having the initial 
value: 

^{To)=Mj^f with j<N, (149) 

and J\fj := ('?/']°\ '?/']°^)+^^^. For these choices of the initial state vector, Cn^ = ioi n > N and 
(|147|) simplifies as 

TV 

V;(r)~5^c„e-^(^-^»)'^"^„ (150) 



n=l 
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Furthermore, for Z = 0, iP{tq) corresponds to a stationary state with a vanishing position 
expectation value for all r G M. For Z > 0, ip{To) does not represent a stationary state and the 
position expectation value is a nonconstant function of time. To determine this function, we 
use the position representation of the state. 

Let \E'(x; r) := (■C^^'') (''"))+ = {'^\p\'4'{t)) be the position wave function for the state vector 
?/'(r). Then, in light of (jHl, (IH3), (EZD, and (ITHn|l . we have 



N 

N 



vl/(x; r) ^ 5^ c„ e-^(---°)'^" r„(x), (151) 



n=l 

where for all ri < 

N 

r„(x) := (x|p|^„) ~ + 



N s-^ ... ^1 /9 ,. , , , . "7rg 

q,k=l 



fl52) 



Next, we employ p4()|l to compute the position expectation value for iP{t). In view of (^SJ), 

N 



(V^(r),X^(r)). 



/ x|^(x;r)|2rfx^ J2 0mne-'("-"«)(''"-^'"\ (153) 

~^ m.n=l 



where 



fi 



1 

N N 



fc,g,D=l "=i 



Figs. El and El show the trajectories traced by the position expectation value ()153|1 in time, 
for the initial state vector ()149|1 with j = 1,2, ■ ■ ■ ,7, Z = 1 and = 10. 



6 Discussion and Conclusion 

In this article we have outlined a general formulation of PT-symmetric (and more generally 
pseudo-Hermitian) quantum mechanics paying attention to the physical aspects of the theory. 
This formulation is consistent with the requirements of quantum measurement theory and al- 
lows for the determination of the physical observables. In fact, to the best of our knowledge, 
this paper is the first to offer an explicit calculation of observables and concrete physical quan- 
tities for a PT-symmetric system with an infinite-dimensional Hilbert space. Perhaps more 
importantly, it proposes a method to identify an underlying classical Hamiltonian that satisfies 
the usual postulates of classical mechanics and a quantization scheme that relates the latter 



37 



n=l n=2 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

X X 



Figure 12: Graph of the trajectory traced by the position expectation value in time for 
the initial state vector tp{0) = Mjip^p with j = 1,2, where Z = 1, = 10, tq = 0, 
T G [0, IG/tt]. The horizontal and vertical axes respectively represent (?/'(r), X'?/'(r))+ and 
r. The r-axis also corresponds to the trajectory for the Hermitian case {Z = 0). Notice 
that IG/tt ^ 5.1 is twice the characteristic period 27c/Ef'^ for the ground state of the 
corresponding Hermitian square well. 
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Figure 13: Graph of the trajectory traced by ('?/'(r), X^/'(r))+ for ^p{0) = -Afj^j with 
j = 3, 4, ■ ■ ■ , 7 and the same parameters and conventions as in Fig. ^1 Note that the 
range of values of the horizontal axis is reduced to amplify the behavior of the trajectories. 
The envelops seem to have the same period as that of Ef'\ i.e., 27r/E|^°^ ^ S/tt. 
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to the defining Hamiltonian of tlie tlieory. We view tliis as a necessary step toward a clearer 
understanding of the potential physical applications of PT-symmetric quantum mechanics. 

Another important outcome of our investigation is that we are now able to consider the 
addition of the interaction terms to a PT-symmetric Hamiltonian H without disturbing the 
structure of its Hilbert space. This is simply done by selecting the additional interaction terms 
from among physical observables. 

Our investigation of the PT-symmetric square well revealed the fact that the underlying 
classical Hamiltonian for this system coincides with that of the ordinary Hermitian infinite 
square well. In other words, the non-Hermiticity effects are quantum mechanical in nature. This 
can be traced back to the simple observation that the non-Hermiticity of the Hamiltonian ()52j) 
only affects the low-lying energy levels. ^'^ 

Our general results confirm the assertion that as a fundamental theory PT-symmetric quan- 
tum mechanics is both mathematically and physically equivalent to conventional quantum me- 
chanics ^Tj. In fact it is this very equivalence that allows for the computation of the physical 
observables. This in turn leads to the natural question whether there is any valid motivation for 
further development of PT-symmetric and pseudo-Hermitian quantum mechanics. Our answer 
to this question is in the affirmative. It is supported by the following observations. 

1. As we showed for any PT-symmetric Hamiltonian if, there is a corresponding Hermi- 
tian Hamiltonian. But the latter is a generically nonlocal (pseudo-differential) operator. 
Therefore, if one is interested in calculating physical quantities that make explicit use 
of the Hamiltonian one is naturally inclined to make use of the original Hamiltonian H 
and the inner product (■, ■)+. However, if one wishes to compute quantities involving the 
position and momentum of the system, then one is essentially forced to use the Hermitian 
picture. In summary, developing PT-symmetric quantum mechanics opens up the possi- 
bility of treating quantum systems with certain nonlocal Hermitian Hamiltonians. The 
classical Hamiltonian for such a system is a real analytic function of x and p that involves 
arbitrarily high powers of p. The pseudo-Hermitian quantization scheme introduced in 
this paper provides a description of the quantum systems associated with these compli- 
cated Hamiltonians. It yields a PT-symmetric quantum Hamiltonian operator that is a 
local (differential) operator. In a sense the use of the PT-symmetric quantum mechan- 
ics is equivalent to trading a complicated nonlocal Hermitian Hamiltonian with a local 
PT-symmetric Hamiltonian.^^ 

^°We do not claim that this is a common feature of all the PT-symmetric quantum Hamiltonians. In general, 

we expect the non-Hermiticity of the quantum Hamiltonian to affect the underlying classical Hamiltonian. We 

leave a more detailed study of this issue for a future publication. 

^^The characterization of all the nonlocal Hermitian Hamiltonians that may be mapped to PT-symmctric or 
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2. The basic ideas so far developed within the framework of pseudo-Hermitian quantum 
mechanics to assess the structure of PT-symmetric quantum mechanics have some re- 
markable applications in relativistic quantum mechanics j^U 122 IHZ] , quantum cosmology 
IHHj . statistical mechanics j^, and magnetohydrodynamics [lUj. This strengthens 
the belief that the study of PT-symmetric quantum mechanics may lead to some con- 
crete advances in other research areas. It is needless to mention the possibility that the 
field theoretical extension of such a study may actually turn out to achieve some of the 
ambitious goals described in [TSj . 

Next, we wish to elucidate the relationship between our approach and the formulation of 
the PT-symmetric quantum mechanics based on the so-called charge-conjugation operator C 
as outlined in ^H]- See also j^. In our approach the metric operator 77+ plays the same role 
as the operator C. In fact as shown in C may be expressed in terms of 77+ according to 

C = 7]l^P. (155) 

Although both 77+ and C determine the inner product of the physical Hilbert space, the expres- 
sion of the latter in terms of 17+ is slightly simpler. It also agrees with the standard mathematical 
approach used in dealing with different inner products on the same vector space. Furthermore, 
the recent attempts |12] at approximate calculations of C for PT-symmetric potentials of the 
form /i^x^ — \^{ix)^ , with /i, A G M, has revealed the remarkable fact that these calculations 
simplify enormously provided that one first computes 77+ and then uses ()155p to determine C.^^ 
This provides a practical justification for the assertion that 77+ is a more basic ingredient of the 
theory than C. The construction of the observables provides a much more concrete evidence 
for the validity of this assertion. Note also that the formulation of the theory that uses C and 
avoids any explicit mention of 77+ leads to the same general conclusions such as the physical 
equivalence of the PT-symmetric and conventional quantum mechanics. A mathematically 
rigorous proof of this statement is given in 

Finally, we wish to comment on whether one can apply the general scheme offered by pseudo- 
Hermitian quantum mechanics to PT-symmetric systems defined on a complex contour. The 
negative attitude expressed by some of the workers on this issue is based on the argument that 
for these systems the eigenfunctions of the Hamiltonian do not belong to the Hilbert space L^(M) 
and hence one cannot define a metric operator 77+ and apply the results of the theory of pseudo- 
Hermitian operators. The problem with this argument is that nowhere in the formulation of 

pseudo-Hermitian Hamiltonians of the standard (kinetic+potcntial) form is a difficult open problem. 

^^What authors of 42 do is to express 77+ as e^'^ ^ calculate Q approximately, and express C as C = e^P. 
Note that this equation is identical with Eq. p55(l that was initially derived in . The calculation of Q makes 
use of the fact that C is a symmetry generator. The observation that operators of the form ?7j~^ry2 generate 
symmetries of a Hamiltonian that is both 771- and 772-pseudo-Hermitian was initially reported in ^j. 
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the pseudo-Hermit ian quantum mechanics ^2] does one assume that the initial Hilbert space 
Ti is L^(M). As explained in Sec. 2, Ti is constructed in two steps: (i) One takes the span Vh of 
the eigenvectors ipn of H (that is assumed to have a real and discrete spectrum) and endows it 
with some arbitrary (positive-definite) inner product, (ii) One performs the Cauchy completion 
of this inner product space to obtain the Hilbert space Ti. The only important condition to be 
checked is whether Ti is separable. This follows from the following simple argument (see also 
[33].) The set {V'n} of the eigenvectors spans Vh, and as Ti is the Cauchy completion of Vh, 
Vh is dense in Ti. This implies that Ti is the closure of the span of {ipn}- Being eigenvectors 
with different eigenvalues, are also linearly independent. Hence {ipn} is a countable basis 
of Ti. In particular, performing Gram-Schmidt orthonormalization on {^„}, one can construct 
a countable orthonormal basis of Ti. This is equivalent to the statement that 7i is a separable 
Hilbert space JHl- This general argument shows that indeed there is no obstruction to employ 
pseudo-Hermitian quantum mechanics to systems defined on a complex contour. The practical 
difference with systems defined on the real axis is that one cannot make a direct use of the 
familiar L^-inner product. It turns out that this does not lead to any insurmountable difficulty 
either. On the contrary, the use of the machinery of pseudo-Hermitian quantum mechanics 
in describing PT-symmetric systems defined on a complex contour has both practical and 
conceptual advantages 
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A Appendix 



In this appendix we present a general method for checking the Hermiticity of the Hamihonian 
h using its power series expansion This provides an interesting explanation for the 

resemblance of the graph of the function $5[5o(x)] to that of the zv(x) as shown in Fig. 
Using p09|) and the fact that v(x) is imaginary, we have 



ht^p2-v(x) + ^p^5,(x)*. (156) 

£=0 

Substituting and in h^ = h, we then find 

^ CO 

v(x) + - - P'^^W*] ^ 0- (157) 

i=Q 

Our purpose is to express the left-hand side of this relation as a power series in p with all 
x-dependent coefficient appearing to the left of powers of p, i.e., obtain a set of functions fe 
such that 

oo ^ oo ^ oo 

J2 /Kx) = v(x) + - - p'Seix)*] = v(x) + - p'U^Y]- (158) 

Then the condition p57|l takes the form 

N 

fe{x) ^ for all £ G {0, 1, 2, ■ ■ ■ }. (159) 

In order to compute fi, we use the following useful identity which may be proven by induc- 
tion on i and use of [x, p] = i. 



k=0 



where / : M — > M is an £-times differentiable function and ^ ^ y k\{£-ky.' 

Using ()160|) . ()158|) . and the fact that v(x) is imaginary, we can express the condition ()159|) . 
after some rather lengthy algebra, as follows. 



For £ = 0, /o ^ yields: 
where 

uo(x) 
wo{x) 



N N 

uo(x) ^ 0, wq(x) ^ iv(x), (16V 



°° Jk 



dx'' 

k=l 



^So{x)] + ]-J2(^-l)'^^hix)]. (163) 



2^' ' c/x^ 

k=l 
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For odd values of £, ~ yields: 



iV 

(x)], we-{x)^0, for all ^ G {1,3,5,- ■■}, 



where 



k=l 

oo 



k=l 



ji + k) 
i\ k\ 



— ^Se+kix)]. 



where 



oo 



k=l 
^ oo 



ii + k)\ 

a k\ 



k=l 



\l + k)\ 
i\ k\ 



— '^[5t+k{x)]. 



(164) 



(165) 
(166) 



For positive even values of £, ~ yields: 

N N ~ 

u,+ (x)^0, w,+ {x) ^ '^[5t{x)l for all £g {2,4,5,- ■■}, (167) 



(168) 
(169) 



Unfortunately, the infinite series appearing in the above relations involve arbitrarily high 
order derivatives of 6i. This reduces their convergence rate appreciably (compared to the series 
expansion ()109|) for h) and amplifies the approximation errors considerably, thus rendering a 
numerical verification of ()16H1 . ()167jl . and ()164j) intractable. However, the second condition in 
(jl61|) provides an interesting explanation for the particular shape of 55 [5o]- Neglecting all the 
terms involving with £ > 0, this condition reads. 



55[5o(x)] ^ ?v(x) 



-Z for -1 < X < 
Z for < X < 1 



This is in remarkable good agreement with graph of 55[5o(x)] as depicted in Fig. El 
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